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We study the interaction of an electromagnetic field with a non-absorbing or absorbing dispersive 
sphere in the framework of complex angular momentum techniques. We assume that the dielectric 
function of the sphere presents a Drude-like behavior or an ionic crystal behavior modelling metallic 
and semiconducting materials. We more particularly emphasize and interpret the modifications 
induced in the resonance spectrum by absorption. We prove that "resonant surface polariton modes" 
are generated by a unique surface wave, i.e., a surface (plasmon or phonon) polariton, propagating 
close to the sphere surface. This surface polariton corresponds to a particular Regge pole of the 
electric part (TM) of the S matrix of the sphere. From the associated Regge trajectory we can 
construct semiclassically the spectrum of the complex frequencies of the resonant surface polariton 
modes which can be considered as Breit-Wigner-type resonances. Furthermore, by taking into 
account the Stokes phenomenon, we derive an asymptotic expression for the position in the complex 
angular momentum plane of the surface polariton Regge pole. We then describe semiclassically the 
surface polariton and provide analytical expressions for its dispersion relation and its damping in 
the non-absorbing and absorbing cases. In these analytic expressions, we more particularly exhibit 
well-isolated terms directly linked to absorption. Finally, we explain why the photon-sphere system 
can be considered as an artificial atom (a "plasmonic atom" or "phononic atom") and we briefly 
discuss the implication of our results in the context of the Casimir effect. 



PACS numbers: 78.20.-e, 41.20.Jb, 73.20.Mf, 42.25.Fx 
I. INTRODUCTION 

In a recent article^, we introduced the complex angular 
momentum (CAM) method in the context of scattering of 
electromagnetic waves from non-absorbing metallic and 
semiconducting cylinders. This allowed us to completely 
describe the resonant aspects of the problem in the fre- 
quency range where the dielectric function is negative. 
In the present article, we shall extend this analysis to 
non-absorbing as well as absorbing metallic and semi- 
conducting spheres. For non-absorbing media, the tran- 
sition from two dimensions to three dimensions induces 
some additional technical difficulties as well as the ap- 
pearance of curvature correction terms in some analytic 
formulas describing surface polaritons (SP's) while the 
physical interpretations remain quasi-identical in form. 
As a consequence, we will pass quickly on certain aspects 
lengthily analyzed in Ref. 1. By contrast, for absorbing 
media, different physical phenomenons occur. They can 
be interpreted from new correction terms (mainly due to 
the imaginary part of the complex dielectric constant) in 
the analytical expression describing the damping of SP's. 

The theory of the resonant surface polariton modes 
(RSPM's) supported by metallic and semiconducting 
spheres has been already studied in numerous works (see, 
for example, Refs. 12.3 4,5.6.7,8.9i ll0lllllll2l and references 
therein) and is currently the subject of a renewed in- 
terest in the fields of nanotechnologies and plasmonics 



(for a recent comprehensive review on this subject, we 
refer to Ref. [H). In the present article, by using CAM 
technique a 14 ! 15 ' 16 in connection with modern aspects of 
asymptotic a 17 ' 18 ! 19 ' 20 , we shall look further into this sub- 
ject. The CAM method has been extensively used in 
physics (we refer to the Introduction of Ref. |l| for a de- 
scription of this method and for a short bibliography). 
We have recently introduced it in the context of electro- 
magnetism of dispersive media (see Refs. [H I21II22I ) . Here, 
by using this method, we shall provide a clear physical 
explanation for the excitation mechanism of the RSPM's 
of the sphere as well as a simple mathematical description 
of the unique surface wave, i.e. of the so-called surface 
polariton (SP), that generates them. 

Our paper is organized as follows. Section II is de- 
voted to the exact theory: we introduce our notations, 
we provide the expression of the S matrix of the system 
and we then discuss the resonant aspects of the prob- 
lem. In Sec. Ill, by using CAM techniques, we establish 
the connection between the SP propagating close to the 
surface sphere and the associated RSPM's. In Sec. IV, 
we describe semiclassically the SP by providing analytic 
expressions for its dispersion relation and its damping. 
Finally, in Sec. V, we conclude our article by explaining 
why the photon-sphere system can be viewed as an ar- 
tificial atom (a "plasmonic atom" or "phononic atom") 
and by briefly discussing the implication of our results in 
the context of the Casimir effect. 
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II. EXACT S MATRIX AND SCATTERING 
RESONANCES 

A. General theory 

From now on, we consider the interaction of an electro- 
magnetic field with a metallic or semiconducting sphere 
with radius a which is embedded in a host medium of 
infinite extent. In the usual spherical coordinate system 
(r, 9, ip) the sphere occupies a region corresponding to the 
range < r < a (region II) while the host medium cor- 
responds to the range r > a (region I). In the following, 
we implicitly assume a time dependence in exp(— iuit) 
for the electromagnetic field and we denote by e c (ui) the 
frequency-dependent dielectric function of the sphere and 
by e/j the constant dielectric function of the host medium. 
Furthermore, we shall use the wave numbers 

k I (uj)=(^j^ and fc n H=(^)^R (1) 

in order to describe wave propagation in regions I and II 
(here c is the velocity of light in vacuum). 

As far as the dielectric function of the sphere is con- 
cerned, we assume it presents a Drude-like behavio r 23 ' 24 



1 



UJ 1 + ijLUpUJ 



or an ionic crystal behavior 2 ' 23 ' 24 



uj\ — uj 2 — i^jUJTU \ 



UJ — IJUJTUJ 



(2) 



(3) 



In both cases, is the high-frequency limit of the dielec- 
tric function and 7 is a phcnomcnological damping factor. 
In Eq. @, Up is the plasma frequency. In Eq. J3]), ujt 
and iv l respectively denote the transverse-optical-phonon 
frequency and the longitudinal-optical-phonon frequency. 
In the first case, SP's follow from the coupling of the 
electromagnetic wave with the plasma wave and are usu- 
ally called surface plasmon polaritons. In the second one, 
SP's follow from the coupling of the electromagnetic wave 
with the longitudinal and transverse acoustic waves and 
are usually called surface phonon polaritons. Eq. ([2]) can 
be used to describe the dielectric behavior of certain met- 
als and semiconductors (Si, Ge, InSb, . . . ) while Eq. ([3]) 
can be used to investigate the optical properties of other 
semiconductors such as GaAs. 

In the following, we shall often consider separately the 
real and imaginary parts of the dielectric function. We 
can write 



with 



e c (w) = + m£M 
e' c M = £o 

4'M = to 



iu 2 p {ujp/u}) 

UJ 2 + (jLOp) 2 



(4) 

(5a) 
(5b) 



(uj 2 - uj 2 )(uj 2 



for the Drude-like behavior and 



Id 2 ) + ("fOJTljj)' 2 



^ujtuj(uj 2 l — uj t ) 



2 ) 2 + (~fUJ T Uj) 2 

t) 



(uj\ - LU 2 ) 2 + ^UJTUjf 



(6a) 
(6b) 



for the ionic crystal behavior. It is important to note that 
the phenomenological damping factor 7 is always smaller 
than the other parameters involved in the expression of 
e c (ui). As a consequence, we can always consider that 

K'MI « 

The S matrix of the sphere is of fundamental impor- 
tance because it contains all the information about the 
interaction of the electromagnetic field with the sphere. 
It can be obtained from Maxwell's equations and usual 
continuity conditions for the electric and magnetic fields 
at the interface between regions I and I I 15 ' 16 i 25 . For our 
problem, the elements of the electric part (TM) of the S 
matrix are given by 



Sf(w) = 1 - 2af(u) 



with 



af(tv) 



(7) 



(8) 



where Cf(ui) and Df(u) are two 2x2 determinants 
which are explicitly given by 

Cf (w) = fc n (w) ipz [k 11 (lu) a] ip'i [k 1 (uj) a] 

-k 1 (ui) ip t [k 1 (w) a] ip' t [k n (w) a] (9a) 

Df{uj) = k u (w) ^ [k n (uj) a] <f >' [jfci ( w ) ] 

-k 1 (uj) C ( e 1] [k 1 (lu) a] i>' t [A: 11 (uj) a] (9b) 

while the elements of its magnetic part (TE) are given 
by 



with 



Sf(u) = 1 - 2af (a,) 
Cf(u) 



af(uj) 



Df(uj) 



(10) 



(11) 



where Cg(uj) and Df 1 (uj) are also two 2x2 determinants 
which are now explicitly given by 

Cf(uj) = k 11 (uj) ip t [k l (uj) a] ip' t [fc n (uj) a] 

-k 1 (uj) i/j e [k 11 (uj) a] if/ t [k 1 (uj) a] (12a) 
Df (uj) = k 11 (uj) ( { e 1] [k 1 (uj) a] i>\ [fc n (uj) a] 

-k 1 (uj) iJt [k u (uj) a] Ci 13 ' [fc 1 (uj) a] . 

(12b) 

In Eqs. © and (p~5|) . we use the Ricatti-Bessel func- 
tions i\ji(z) and (z) which are linked to the spherical 
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Bessel functions je(z) and h!f\z) by ipe(z) = zji(z) and 

C,f\z) = zh ( e 1] (z) (see Ref.EjJ. 

From the S matrix elements, we can in particular con- 
struct the scattering cross section er sca and the absorp- 
tion cross section cr a b s of the spher o 15 i 16 . They can be 
expressed in terms of the coefficients af(u>) and af^(id) 
and they are given by 



Csca(w) 



2tt 



£(2^+l)[|afM| 2 + |af(o,)| 2 ] (13) 



and 



Cabs(w) 



2tt 



[k l (u)] J £J 



^(2£ + l){[Reaf(id)-\af(id)\' 



+ [ReafM-|afH| 2 ]}.(14) 



From the S matrix elements, we can also precisely de- 
scribe the resonant behavior of the sphere as well as the 
geometrical and diffractive aspects of the scattering pro- 
cess. Here the dual structure of the S matrix plays a 
crucial role. Indeed, the S matrix is a function of both 
the frequency id and the angular momentum index I. It 
can be analytically extended into the complex w-plane as 
well as into the complex A-plane or CAM plane. (Here A 
denotes the complex angular momentum index replacing 
£+1/2 where £ is the ordinary momentum inde x 1 W 6 .) 
The poles of the S matrix lying in the fourth quadrant 
of the complex w-plane are the complex frequencies of 
the resonant modes. These resonances are determined 
by solving 



D 



E,M 



(w)=0 for £=1,2,3,. 



The solutions of ([T5|l are denoted by oj£ p 



,(°) 



(15) 

iTi p /2 

where u>^ > and T( p > 0, the index p permitting us to 
distinguish between the different roots of (|15[) for a given 
L In the immediate neighborhood of the resonance u>£ p , 
(ui) has the Breit-Wigner form, i.e., is proportional 



to 



E,M 



,(0) 



iTip/2 



(16) 



The structure of the S matrix in the complex A-plane 
allows us, by using integration contour deformations, 
Cauchy's Theorem and asymptotic analysis, to provide 
a semiclassical description of scattering in terms of rays 
(geometrical and diffracted). In that context, the poles 
of the S'-matrix lying in the CAM plane (the so-called 
Regge poles) are associated with diffraction. They are 
determined by solving 



r\E,M I \ 



for id > 0. 



(17) 



Of course, when a connection between these two faces of 
the S matrix can be established, resonance aspects are 
then semiclassically interpreted. 



In the following, we shall present some numerical re- 
sults. We have chosen to restrict ourselves to particular 
configurations, i.e., to particular values of the parameters 
eh, £oo, id p , lutci/c, (dLa/c and 7. These values are phys- 
ically realistic or, more precisely, of the same order than 
physically realistic values (see, for example, Refs. I2II23T). 
In fact, even for different configurations, the results we 
have numerically obtained and that we shall discuss in 
this paper remain valid. 



B. Sphere with a Drude-like behavior 

In Figs. [1] and O we consider the resonant aspects of 
a sphere embedded in vacuum (e^ = 1) and we assume 
that its dielectric function presents the Drude-like behav- 
ior given by Eq. ^ with = 1 and id p a/c = 2ir. Wc 
examine both the non-absorbing case with 7 = in Fig. [1] 
and the absorbing case with 7 = 1/100 in Fig. [2] In the 
non-absorbing case, we display the scattering cross sec- 
tion in Fig. [TJl and, in the absorbing case, we display the 
absorption cross section in Fig. [2ji. These cross sections 
are both plotted as functions of the reduced frequency 
tda/c. On the two figures, rapid variations of sharp char- 
acteristic shapes can be observed. This strongly fluctuat- 
ing behavior is due to scattering resonances: when a pole 
of the S matrix is sufficiently close to the real axis in the 
complex cj-plane, it has a strong influence on the cross 
section [see Eq. (fT6"|) ]. In Figs.Qi) and[^b, resonances are 
exhibited for the two configurations previously consid- 
ered. A one-to-one correspondence between the peaks of 
the cross sections and the resonances near the real ua/c- 
axis can be clearly observed in certain frequency ranges. 

More precisely and more generally, for the dielectric 
function @ there exists in the frequency range where 
e' c (td) < (i.e., where lj < lo p if we neglect terms in 
7 2 ) a family of resonances associated with S E (TM res- 
onances). They are close to the real axis of the complex 
w-plane and they converge, for large £, to the limiting 
complex frequency io s satisfying 



,(w s ) + e h = 0. 



(18) 



The real and imaginary parts of td s are easily found per- 
turbatively by inserting 

id s = id' s + iid" (19) 

into (HHJ) and by taking into account ([3]) with |e"(a>)| <C 
|e£(w)|. By assuming that |co>"| <C id' s and by using a first- 
order Taylor series expansion of e c (u> s ), we find that id' s 
must satisfy 



and that 



e' c (id' s ) + e h = 0. 



d Re e' c (id) / did 



(20a) 



(20b) 
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FIG. 1: a) Scattering cross section a BC a. b) Scattering res- 
onances in the complex wa/c-plane. We consider a non- 
absorbing sphere: e c (ui) has the Drude type behavior with 
eoo — 1, uj p a/c — 2tt and 7 = while en = 1. Dots (■) corre- 
spond to poles of S E {to) while plus (+) correspond to poles 
oiS M {tu). 



FIG. 2: a) Absorption cross section <7 a b s - b) Scattering reso- 
nances in the complex wa/c-plane. We consider an absorbing 
sphere: e c (uj) has the Drude type behavior with — 1, 
uj p a/c = 2ir and 7 = 1/100 while eh = 1. Dots (■) corre- 
spond to poles of S e (lu) while plus (+) correspond to poles 
of S m (lo). 



By using Eqs. (j5a|) and j5b|) and by neglecting terms in 
7 2 , we then obtain 



^/l + eh/eoo 



The general formulas (|21a[) and (|21b|) describe very 
well the accumulation of resonances observed in Figs.[Tj3 
and [2b in the frequency range where e' c (u>) < 0. It is 
important to note the existence of a shift in the imagi- 
nary part of the resonance spectrum [see Eq. (|21b[) ]. It is 
well highlighted in Fig. [2Jd. It is associated with absorp- 
tion and proportional to the phenomcnological damping 
factor 7. At first sight, it can appear surprising that it 
does not depend on the dielectric constants and eh- 
In fact, such a dependence only appears by working at 
higher orders in the perturbative expansion used. 

The resonances we have just considered are associated 
with the electric part S E of the S matrix and they corre- 
spond to TM resonant modes of the sphere whose excita- 
tion frequencies belong to the frequency range in which 
e' c (u>) < 0. From now on (i.e., in Sees. Ill and IV), we 
shall more particularly focus our attention on the phys- 
ical interpretation of these resonant modes. We shall 
prove that the corresponding resonances are generated 
by an exponentially attenuated SP propagating close to 
the sphere surface, this fact justifying the term RSPM's 
used to denote the associated resonant modes. Further- 
more, we shall show that the SP is the analog, in the 
large radius limit, to that supported by the plane inter- 
face. Of course, there also exists, in the whole frequency 



range, other families of resonances with "higher" imagi- 
nary parts and which are associated with both the elec- 
tric and the magnetic parts of the S matrix. They corre- 
spond to the excitation of TM and TE resonant modes of 
the sphere. From a physical point of view, these modes 
are less interesting due to their shorter lifetime. Indeed, 
in the new field of plasmonics, those are the SP's with 
long propagation lengths and therefore very small atten- 
uations that are especially interesting from the point of 
view of practical applications. Furthermore, if we con- 
sider the system photon-sphere as an artificial atom for 
which the photon plays the usual role of the electron (a 
point of view we shall push farther in Sees. Ill and IV), 
we must then keep in mind that, in the scattering of a 
photon with frequency u>^ , a decaying state (i.e., a qua- 
sibound state) of the photon-sphere system is formed. It 
has a finite lifetime proportional to l/Ti p . The resonant 
states whose complex frequencies belong to the family 
generated by the SP are therefore the most interesting 
because they are very long-lived states. 



C. Sphere with an ionic crystal behavior 

In Figs. [3] and [4l we consider the resonant aspects of a 
sphere embedded in vacuum (e^ = 1) and we assume that 
its dielectric function presents the ionic crystal behavior 
given by Eq. ^ with = 2, lutcl/c — 2tt and Wi,o/c = 
3ir. We examine both the non-absorbing case with 7 = 
in Fig.[3]and the absorbing case with 7 = 1/100 in Fig. 2] 
In the non-absorbing case, we display the scattering cross 
section in Fig. [3^ and, in the absorbing case, we display 
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coa/c 



FIG. 3: a) Scattering cross section a BC a. b) Scattering res- 
onances in the complex wa/c-plane. We consider a non- 
absorbing sphere: e c (fJ) has the ionic crystal behavior with 
too = 2, (jJtci/c — 2n, lolcl/c = 37r and 7 = while eh = 1. 
Dots (•) correspond to poles of S (u>) while plus (+) corre- 
spond to poles of S M (cj). 




FIG. 4: a) Absorption cross section <7 a b s . b) Scattering reso- 
nances in the complex wa/c-plane. We consider an absorbing 
sphere: e c (u>) has the ionic crystal behavior with e^ = 2, 
lotcl/c — 2-7T, ujl<i/c = 3-7T and 7 = 1/100 while eh = 1. Dots 
(•) correspond to poles of S E (ui) while plus (+) correspond to 
poles of S M (co). 



the absorption cross section in Fig. 0}t. In Figs. [5Jd and 
2Jd, resonances are exhibited for the two configurations. 
A one-to-one correspondence between the peaks of the 
cross sections and the resonances near the real wa/oaxis 
can be clearly observed in certain frequency ranges. 

Here again, we focus our attention on the resonances 
existing in the frequency range where e' c (u)) < (i.e., 
where lot ^$ w P$ if we neglect terms in j 2 ). There 
is a family of resonances associated with S E (TM res- 
onances) close to the real axis of the complex cj-plane. 
They converge, for large i, to the limiting complex fre- 
quency lo s still satisfying Eq. (|18[) with e c now given by 
Eq. ^ . The real and imaginary parts of oj s = lu' s + ico'J 
can be obtained perturbatively and we have 



'cj| + (e/i/eoo)^ 



1 + eft/e 



7Wy 



(22a) 
(22b) 



The general formulas (|22ap and (|22b|) describe very 
well the accumulation of resonances observed in Figs. [3)3 
and [Dd in the frequency range where e' c (u>) < as well as 
the shift in the imaginary part of the resonance spectrum. 

In Sees. Ill and IV, we shall provide a physical inter- 
pretation of the family of resonances we have just con- 
sidered and we shall prove that these resonances are gen- 
erated by an exponentially attenuated SP propagating 
close to the sphere surface and analog, in the large radius 
limit, to that supported by the plane interface. Of course, 
there also exists other families of resonances but we shall 
not focus our attention on them even if they seem to us 
physically interesting due to their rather "low" imaginary 



parts (see the frequency range wa/c < ujTa/c in Figs. [SJd 
andUJa). In fact, it is not possible to interpret them in 
term of SP's analog to that supported by the plane in- 
terface because, in frequency ranges where e' c (w) > 0, 
no surface wave can be supported by the plane interface. 
We think however that a semiclassical description could 
be achieved (maybe in terms of whispering-gallery-type 
surface waves) but it is out of the scope of the present 
work. 



D. Comparison between cylinders and spheres 

We conclude this section by making a brief comparison 
with the results obtained in our previous study concern- 
ing metallic and semiconducting cylinders^. Of course, 
in Ref. [l| we have only considered non-absorbing cylin- 
ders. As a consequence, a comparison between cylin- 
ders and spheres can be achieved only in this restricted 
context and we can then notice that the cross sections 
for the cylinders and the spheres as well as the spectra 
of resonances are rather similar. However, it should be 
noted that in the scattering by a sphere both the TM 
and TE polarizations contribute to the cross section (we 
recall that for the cylinder, the two polarizations can be 
studied separately) and it is important to note that, for 
both scatterers, SP's and their associated RSPM's corre- 
spond to only one polarization (the TE polarization for 
the cylinder and the TM polarization for the sphere) . Of 
course, these last results remain valid even in the pres- 
ence of absorption. 



6 



III. SEMICLASSICAL ANALYSIS: FROM THE 
SP REGGE POLE TO THE COMPLEX 
FREQUENCIES OF RSPM'S 

As already mentioned in Sec. II (see also Refs. Hil . 
ITEllTil ). in the CAM approach, Regge poles determined 
by solving Eq. (JT7J) are crucial to describe diffraction 
as well as resonance phenomenons in terms of surface 
waves. From the Regge trajectory associated with the 
SP supported by the metallic or semiconducting sphere, 
i.e., from the curve Asp = Asp(w) traced out in the CAM 
plane by the corresponding Regge pole as a function of 
the frequency, we can more particularly deduce: 
(i) the dispersion relation 



A;sp(w) 



Re Asp(w) 



(23) 



of the SP which connects its wave number kgp with the 
frequency uj, 

(ii) the damping ImAsp(^) of the SP, 

(iii) the phase velocity v p as well as the group veloc- 
ity v g of the SP given by 

uia d toa , 

Wp = M^) ^ Vg ^ d Re A S p M ' (24) 

(iv) the semiclassical formula (a Bohr-Sommerfeld 
type quantization condition) which provides the location 

of the excitation frequencies wjgp of the resonances gen- 
erated by the SP: 



ReA SP (4sp) =* + l/2 €=1,2,. 



(25) 



(v) the semiclassical formula which provides the 
widths of these resonances 



£SP 



Im Asp (lo) (e?Re Asp (w) / doj) 



(dReX SP (uj)/duj) 2 + (dim. \ S p(uj)/dtu) 



and which reduces to 

r^sp 



Im Asp(w) 



d Re \sp(uj)/du; 



(26) 



(27) 



in the frequency range where the condition 
\dRe \sp(uj)/duj\ » \dlm Xsp(ui)/dto\ is satisfied. All 
these results can be established by generalizing, mutatis 
mutandis, our approach and our calculations developed 
in Refs[l|[2l| for dispersive cylinders. The transition 
from the dimension 2 to the dimension 3 induces some 
additional technical difficulties (vectorial treatment, ex- 
istence of a caustic, asymptotics for spherical harmonics 
...) which can be rather easily overcome following and 
extending the works of Newton in quantum mechanics 
(see Ch. 13 of Ref. Il4h and the works of Nussenzveigi^ 
and Grandyi^ in electromagnetism of ordinary dielectric 
media. 

Figs. [5] and [5] exhibit the distribution of Regge poles 
for a sphere embedded in vacuum when e' c (u>) < 0. We 
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FIG. 5: Regge poles in the complex angular momentum plane. 
e c (w) has the Drude type behavior with — 1 and u v a/c = 
27r while th = 1. The distribution corresponds to uia/c — 4 
and we have e' c (<+>) < 0. Dots (■) and crosses (x) correspond 
respectively to Regge poles for 7 = and for 7 = 1/100. 
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FIG. 6: Regge poles in the complex angular momentum plane. 
e c (w) has the ionic crystal behavior with eoo = 2, ura/c = 2ir 
and LOha/c — 3ir while eh = 1. The distribution corresponds 
to uua/c = 8.3 and we have e' c (u>) < 0. Dots (■) and crosses 
(x) correspond respectively to Regge poles for 7 = and for 
7 = 1/100. 



only consider the Regge poles of the electric part S E of 
the S matrix (TM polarization) for the configurations 
numerically studied in Sec. II. These Regge pole distri- 
butions are rather similar to the distributions associated 
with the ordinary dielectric spher o 15 ! 16 . However, some- 
thing new occurs: there exists a well-identified particular 
Regge pole lying in the first quadrant of the A-plane and 
close to the real axis. This new Regge pole Asp(w) is 
associated with the SP orbiting around the metallic or 
semiconducting sphere. Absorption induces an impor- 
tant modification of its imaginary part while it leaves 
unchanged the position of the other Regge poles. For 
other configurations (i.e., for other values of the param- 
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FIG. 7: Regge trajectory for the SP Regge pole: comparison 
between a non-absorbing and an absorbing sphere. e c {u)) has 
the Drude type behavior with = 1 and ui p a/c = 2n while 
£h = 1. As Loa/c — > iu' B a/c, the real part of the SP Regge pole 
always increases indefinitely while its imaginary part vanishes 
for 7 = and increases indefinitely for 7^0. 



-Y=0 
-Y*0 



TM polarization 



o 



6 "v* 7 
coa/c 



Y=0 






_ 7*0 













coa/c 



FIG. 8: Regge trajectory for the SP Regge pole: comparison 
between a non-absorbing and an absorbing sphere. e c {u)) has 
the ionic crystal behavior with eoo = 2, ura/c = 2n and 
ujlcl/c = 3-7T while = 1. As uja/c — > uo' 3 a/c, the real part of 
the SP Regge pole increases indefinitely while its imaginary 
part vanishes for 7 = and increases indefinitely for 7 ^ 0. 



eters e^, e/j, u p , lot, and 7), the Regge pole distri- 
butions are not globally different from those of Figs. 
and[H The SP Regge pole Asp(w) is still present. By 
contrast, when e' c (oj) > or when we consider the TE 
polarization, the SP Regge pole never exists. 

In Figs. [7J and [5J we have displayed the Regge tra- 
jectories of the SP for the two configurations previously 
studied. Absorption does not modify the global behav- 
ior of the real part of the SP Regge pole. It should be 
also noted that, as to — > lj' s , this real part increases in- 
definitely. By contrast, absorption increases significantly 
the imaginary part of the Regge trajectory of the SP: for 



TABLE I: The first complex frequencies of RSPM's (TM po- 
larization). e c (ui) has the Drude type behavior with e^ = 1, 
u> v a/c = 2-7T and 7 = while eh = 1. 





Exact 


Exact 


Semiclassical 


Semiclassical 


I 


(w)gp)o/c 


(r eSP /2)a/c 


(«isp)°/ c 


(r ffiP /2)o/c 


1 


0.897849 


0.407409 


0.944599 


0.381875 


2 


1.762396 


0.391747 


1.734451 


0.380517 


3 


2.475717 


0.268405 


2.430280 


0.267887 


4 


3.019638 


0.136587 


2.992907 


0.138470 


5 


3.412422 


0.049796 


3.404834 


0.050338 


6 


3.683304 


0.012179 


3.682351 


0.012207 


7 


3.862800 


0.001935 


3.862752 


0.001934 


8 


3.982206 


0.000209 


3.982204 


0.000208 


9 


4.065198 


0.000016 


4.065196 


0.000016 


10 


4.125633 


0.000001 


4.125631 


0.000001 



TABLE II: The first complex frequencies of RSPM's (TM 
polarization). e c (uj) has the Drude type behavior with = 
1, LJ p a/c = 2n and 7 = 1/100 while eh = 1. 





Exact 


Exact 


Semiclassical 


Semiclassical 


I 


(u& P )a/c 


(r lSP /2)a/c 


(«isp)a/c 


(r ffi p/2)a/c 


1 


0.894546 


0.408326 


0.942572 


0.382753 


2 


1.757428 


0.395822 


1.731152 


0.384335 


3 


2.470440 


0.276889 


2.426006 


0.276611 


4 


3.015603 


0.149847 


2.988477 


0.152411 


5 


3.410258 


0.067443 


3.401552 


0.068144 


6 


3.682529 


0.033581 


3.681107 


0.033585 


7 


3.862589 


0.026173 


3.863243 


0.026087 


8 


3.982111 


0.026321 


3.983622 


0.026089 


9 


4.065111 


0.027345 


4.067129 


0.027396 


10 


4.125541 


0.028160 


4.128020 


0.028289 



the non-absorbing sphere, this imaginary part vanishes as 
lo — > lo' s while, for the absorbing sphere, it first exhibits 
a minimum and then increases indefinitely as u — > lo' s . 
In other words, absorption has a strong influence on the 
damping of the SP but does not modify its dispersion re- 
lation. For other configurations (i.e., for other values of 
the parameters e^, eh, u) p , ujt, ljl and 7), the SP Regge 
pole behavior is very similar. The minimum of the imag- 
inary part of the Regge trajectory always exists (for the 
absorbing sphere) . This feature could be interesting with 
in mind practical applications using real materials. 

Tables HI [Til Hill and IIVI present samples of complex fre- 
quencies of RSPM's for the two configurations previously 
considered. They are calculated from the semiclassical 
formulas and (|26p by using the Regge trajectories 
numerically determined by solving (jTTJ) (see Figs. [7] and 
HJ. We can observe a very good agreement between the 
"exact" and the semiclassical spectra for "high" frequen- 
cies as well as a rather good agreement for "low" frequen- 
cies. Furthermore, from the behavior of Regge trajecto- 
ries near the limiting frequencies lu' s and the semiclassical 
formulas and Prj)) . we easily obtain the existence of 
the families of resonances close to the real axis of the com- 
plex w-plane which converge, for large t, to the limiting 
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TABLE III: Some complex frequencies of RSPM's (TM po- 
larization). e c (uo) has the ionic crystal behavior with = 2, 
lotcl/c — 2ty, u)lo,/c = 3tv and 7 = while eh = 1. 





Exact 


Exact 


Semiclassical 


Semiclassical 


£ 


(cd e % )a/c 


(r £S p/2)a/c 


(o4sp)a/c 


(r <S p/2)o/c 


9 


7.436965 


0.111396 


7.405909 


0.115670 


10 


7.700893 


0.037894 


7.694243 


0.038461 


11 


7.887534 


0.009632 


7.886710 


0.009713 


12 


8.015334 


0.001839 


8.015239 


0.001844 


13 


8.104140 


0.000274 


8.104087 


0.000273 


14 


8.168509 


0.000034 


8.168477 


0.000034 


15 


8.217128 


0.000005 


8.217069 


0.000003 


16 


8.255038 


0.000002 


8.254956 


0.000000 


17 


8.285324 


0.000001 


8.285231 


0.000000 



TABLE IV: Some complex frequencies of RSPM's (TM po- 
larization). e c {ijj) has the ionic crystal behavior with = 2, 
ujtcl/c = 2-7T, LOLdjc = 37T and 7 = 1/100 while en = 1. 





Exact 


Exact 


Semiclassical 


Semiclassical 


I 


(cd e lp)a/c 


(r ffiP /2)a/c 


(wisp)°/ c 


(L ffi p/2)o/c 


9 


7.433460 


0.131422 


7.399537 


0.137680 


10 


7.699297 


0.060583 


7.690606 


0.061609 


11 


7.887001 


0.034650 


7.885688 


0.035096 


12 


8.015186 


0.028573 


8.016067 


0.028737 


13 


8.104082 


0.028141 


8.105881 


0.028098 


14 


8.168462 


0.028649 


8.170808 


0.027735 


15 


8.217080 


0.029138 


8.219861 


0.029527 


16 


8.254988 


0.029514 


8.258154 


0.029563 


17 


8.285273 


0.029801 


8.288906 


0.029931 



frequency u>' s + iui'J. In conclusion, we have established a 
connection between the complex frequencies of RSPM's 
and a particular surface wave, the so-called SP, described 
by a particular Regge pole of the electric part of the S 
matrix and which orbits around the sphere. 

We conclude this section by making a brief comparison 
with the results obtained in our previous study concern- 
ing metallic and semiconducting cylinders^. Of course, 
such a comparison can be achieved only in the non- 
absorbing case. From Regge trajectories (see Figs. 6 and 
7 of Ref. [0 and Figs. [7J and [5] of the present article), we 
can observe that the behavior of the SP orbiting around 
a metallic/semiconducting sphere is, at first sight, rather 
similar to the behavior of the SP orbiting around a metal- 
lic/semiconducting cylinder, even if they correspond to 
different polarizations (TE polarization for the cylinder 
and TM polarization for the sphere). In fact, as we shall 
see in the next section, in the absence of absorption the 
transition from two dimensions to three dimensions in- 
duces some curvature corrections on the wave numbers 
^Sp(w) of the SP's and the behaviors are identical only in 
the radius limit a — > 00, i.e., in the flat interface limit. In 
the presence of absorption, we shall prove than a supple- 
mentary correction associated with the imaginary part of 
the complex dielectric constant is necessary in order to 
interpret the behavior observed in Figs. [7] and [51 

IV. SEMICLASSICAL ANALYSIS: 
ASYMPTOTICS FOR THE SP AND PHYSICAL 
DESCRIPTION 

An analytical expression for the Regge pole Asp and 
therefore a deeper physical understanding of the SP be- 
havior can be obtained by solving Eq. (jTTJ) for A = 
Asp- For the TM polarization, Eq. (fTTj) reduces to [see 
Eq. j9b}] 

v^H clsp-i/2 (vgWc) _ ^A S p-i/2 (v^H^/c) 

^ Cxl-i/2 (v^<Wc) Va S p-i/2 ( v / ec(w)aw/c) 

(28) 



This equation cannot be solved exactly but only pertur- 
batively. With this aim in view, we must first replace 
Ricatti-Bessel functions by spherical Bessel functions. 
We then use their relations with the ordinary Bessel func- 
tions (see Ref. HI) 



JaO) = J ^J\+i/2{z) and h { ^{z) = ^ ^H^_ 1/2 (z) 

(29) 

and we finally replace the Bessel function J\(z) by the 
modified Bessel function I\(z) (see Ref. l26h in order to 
take into account the fact that Ree c (ti;) < 0. Eq. [)28p 
then reduces to 

1 g&pCVwWc) J_ /_c\ 

= 1 r\ SP W-tc{u)ua/c) 1 /_c_\ 

yj-e c {u) I XsP (y/-e c {u)u)a/c) 2e c {uj) \ujaJ ' 

(30) 

This equation must be compared with Eq. (26) of Ref. 1 
which provides the SP Regge pole for the cylinder. The 
first term on the left-hand side and the right-hand side 
of Eq. ((30)) are exactly those appearing in Eq. (26) of 
Ref. [l|. The two others terms are simple curvature correc- 
tions due to the change of dimension. As a consequence, 
Eq. (|3"0)) can be solved following the method used in or- 
der to solve Eq. (26) of Ref. 0] but, now, it is important 
to take carefully into account the fact that the dielectric 
function e c (ui) has an imaginary part [see Eqs. (|5b[) and 
(|6b[) ]. It should be noted that the existence of such cur- 
vature corrections has been first observed by Berry^i in 
the restricted case of non-dispersive and non-absorbing 
spheres. 

On the right-hand side of (f3"0"]) . by assuming | Asp | 
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\\/— eJuJ)uia/c\, we can writei^ 

1 j a sp (a/- £ c(^)W c ) 
y/-e c (uj) h sp (y/-e c (u;)uja/c) 



[^p- £c MHc) 2 ] 1/2 



6 c (w)(wa/c) 



(31) 



On the left-hand side of (|30|) . by assuming |Agp| ^> 
y/ehUJa/c, we can write ^ 28 ' 29 

1 ff£,'(v^wo/c) 



2 

SP 



e/^wa/c) 2 ] 172 



e h (uja/c) 



x 1 



z e 



(32) 



where 



a(A,z) = (A"-z 2 ) 1/2 -Aln 



A + (A — z 2 ) 



2\l/2 



(33) 



It should be noted that in Eq. (|32|) we have taken into 
account an exponentially small contribution (the term 
exp [2a(Asp, y^wa/c)]) which lies beyond all orders in 
perturbation theory. This term can be captured by care- 
fully taking into account Stokes phenomeno n 17 i 18 ' 19 i 20 
and is necessary to extract the asymptotic expression of 
the imaginary part of Asp . (In Eq. (|32"|) we have given to 
the Stokes multiplier function the value 1/2.) For more 
precision, we refer to our previous article^. 

By using ([31]) and J32|) as well as |e"(w)| < \e' e {u)\, 
Eq. (j3"0|) can be solved perturbatively and we find 



Re 



v c / y e h + e' c {uj) y 



(34a) 



and 

with 

Imi A SP (o;) 



Im Asp(w) = Imx Asp(o;) + Im 2 Asp(w) 



[Re Asp(u)] 2 - th (uja/c) 2 
Re Asp(w) 



l _ € h /e' c (uj) /_c_\ 
2y/-[e h + e' c (u)] W 



Im 2 A S p(w) ~ I — 



ua\ I e h e' c (uj) e h e'^(u;) 



c J y e h + e> c {u) 2e> c (u)[e h + e' c {u)] 



2y/ -[e h + e' c {uj)} Vwfl/ 



(34b) 

exp{2a[Re Asp(w), ^/e^uia/c}}, 

(34c) 
(34d) 



Equations f34afl . (|34b[) . (f34c)) and (|34d[) provide an- 
alytic expressions for the dispersion relation and the 
damping of the SP. The following important features 
must be noted: 

- The SP only exists in the frequency range where 
eh + e' c {u>) < 0. Its dispersion relation [see Eqs. (|34a[> and 
(|23p] only depends on the real part of the dielectric func- 
tion: thus, it is slightly modify by absorption. Its atten- 
uation [see Eqs. J34b|, (|34c[) and (|34d[) ] is compounded 
by two contributions. The first one [see Eq. (|34c| only 
depends on the real part of the dielectric function and 
presents an exponentially small attenuation. It fully de- 
scribes the attenuation of the SP propagating on a non- 
absorbing sphere. The second contribution [see Eq. (|34d| 
is proportional to the imaginary part of the dielectric 
function and is therefore directly linked to absorption. 
It increases indefinitely as u> — > u>' s and semiclassically 



explains the behavior already described in Sec. III. 

- The wave number ksp(u>) associated with the 
SP is obtained from (|34ap and (|23[) and is given by 



«SP 



1 



2y/-[e h + e' c (u})] \uaJ 



(35) 



This relation could permit us to derive analytically the 
phase velocity v p = uj/ksp(u) as well as the group veloc- 
ity v g = dui/dk SP (ui) of the SP [see also Eq. ([24)]. 

- For a — > oo - i.e., in the flat interface limit - we 
recover the usual dispersion relation of a SP supported 
by a flat metal-dielectric or semiconductor-dielectric in- 
terface (see, for example, Ref. [30h - 
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FIG. 9: Regge trajectory for the SP Regge pole. Comparison 
between exact and asymptotic theories. e c (ui) has the Drude 
type behavior with = 1, uj p a/c = 2vr and 7 = while 
e h = l. 



FIG. 11: Regge trajectory for the SP Regge pole. Comparison 
between exact and asymptotic theories. e c (u) has the ionic 
crystal behavior with e x = 2, cj^a/c = 2-7T, cj^a/c = 37r and 
7 = while en ~ 1. 




coa/c 



FIG. 10: Regge trajectory for the SP Regge pole. Comparison 
between exact and asymptotic theories. e c (ui) has the Drude 
type behavior with = 1, ui p a/c = 2n and 7 = 1/100 while 
e h = 1. 



FIG. 12: Regge trajectory for the SP Regge pole. Comparison 
between exact and asymptotic theories. e c (uj) has the ionic 
crystal behavior with = 2, lotcl/c = 2-7T, lulcl/c = 37r and 
7 = 1/100 while e h = 1. 



- The imaginary part (|34cp of Asp vanishes for 
a — > 00: the SP supported by the flat interface between 
an ordinary dielectric and a non-absorbing metallic or 
semiconducting medium has no damping (see, for exam- 
ple, Ref.[I|. 

- By inserting the expression (|34ap for Asp(w) into 
the Bohr-Sommcrfeld quantization condition (|25[) . we ob- 
tain an eighth-order polynomial equation which can be 
solved numerically and which provides rather precisely 
the resonance excitation frequencies w|gp. 

- We have numerically tested the formulas (|34a| . 
(|34b[) . (|34cp and (|34d[) for various values of the parame- 
ters eoo, eh, Up, u>t, Wi and 7. In all cases, they provide 
rather good approximations for Re Asp and Im Asp (see 
Figs. O EB HH and [T3] for the two configurations previ- 
ously studied). 



The connection with the results obtained in our pre- 
vious study concerning metallic and semiconducting 
cylinders 1 can be made. The transition from the cylinder 
to the sphere induces some curvature corrections on the 
real and imaginary parts of the Regge pole of the SP [the 
second term on the right-hand side of (|34a|) . the third 
term on the right-hand side of (|34cji and the fourth term 
on the right-hand side of (|34d|) ] which come from the sec- 
ond terms on the left-hand side and the right-hand side of 
Eq. ([317)) . Without the curvature corrections, Eqs. (|34ap . 
plbj) . (j3"4"c|) and (|34d)) reduce to 




(36a) 
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and 



Im Asp (u>) = Inii Asp (u>) + Ini2 Asp (w) 



(36b) 



with 

Inii A sp (cj) 



[Re A SP (u)] 2 - e h (uja/c) 2 
Re A SP (w) 
x exp{2a[Re Agp(w), y/enuja/c]}, 



Ini' 



SPl >~ v c if £/1 + <H2 £ ' c H[ eft + e j 



(36c) 

Hi' 

(36d) 



These last expressions are the semiclassical formulas pro- 
viding the dispersion relation and the damping of the SP 
propagating on the cylinder in the non-absorbing and 
absorbing cases. They are much more general that the 
corresponding formulas obtained in our previous paper—. 
Indeed, Eq. (|36a[) has been already obtained in Ref. 1 
in the restricted context of the non-absorbing cylinder 
[see Eq. (40a)] and, in the same way, Eq. (|36cp has been 
obtained in a less practical form [see Eq. (40b) in Ref. [l[ . 



V. CONCLUSION AND PERSPECTIVES 

In this article, we have introduced the CAM method 
in order (i) to study the interaction of electromagnetic 
waves with non-absorbing or absorbing metallic and 
semiconducting spheres and (ii) to completely describe 
the resonant aspects of this problem. This allows us to 
provide a physical explanation for the excitation mech- 
anism of RSPM's as well as a simple mathematical de- 
scription of the surface wave (i.e., the SP) that generates 
them. It should be noted that our approach is not lim- 
ited to the metals and semiconductors described by ((2|) 
and © but remains still valid for more general materials 
(see the conclusion of Ref. [l|) . 

In Ref. 0], we have developed a new picture of the 
photon-cylinder system (see also the analysis of the works 
of Ito and Sakod a 31 i 32 in the conclusion of Ref. HI): it 
can be viewed as an artificial atom for which the photon 
plays the role of an electron. Mutatis mutandis, this pic- 
ture is still valid for the photon-sphere system: RSPM's 
are long-lived quasibound states for this atom, the as- 
sociated complex frequencies are Breit-Wigner-type res- 
onances while the trajectory of the SP which generates 



them and which is supported by the sphere surface is 
a Bohr-Sommerfeld-type orbit. Recently, Guzatov and 
Klimov 3 - 3 . have also developed the analogy between a 
metallic sphere and an ordinary atom and, more gen- 
erally, between a cluster of metallic spheres and ordinary 
molecules, introducing on that occasion the terms "plas- 
monic atom" and "plasmonic molecules" . It should be 
noted that their approach uses the quasi-static approx- 
imation for the description of the electromagnetic field. 
As a consequence, they have found that the "energy lev- 
els" of the photon-sphere system are real (bound states). 
In fact, as we have shown, the imaginary parts of the 
"energy levels" do not vanish (quasi-bound states): they 
correspond to exponentially small terms lying beyond all 
orders in perturbation theory. 

Our results could be useful (i) in the context of three- 
dimensional photonic crystal physics (the existence of 
dispersionless band for the arrays of metallic or semicon- 
ducting spheres is associated with the excitation of the 
sphere RPSM's), (ii) in cavity quantum electrodynamics 
and, more generally, (iii) in the context of nanotechnolo- 
gies and plasmonics. Here, we shall more particularly 
focus our discussion to the possible applications to the 
Casimir effect. 

Since the precise experiments carried out by 
Lamoreau?* 3 ^ in 1997, it is necessary to completely de- 
scribe, from a theoretical point of view, the Casimir inter- 
action between two metallic spheres or between a metal- 
lic sphere and a plane (for a review on recent experi- 
ments and problematic, we refer to Ref. 35). A semi- 
classical description of this interaction could be achieved 
by extending to electrodynamics the Korringa-Kohn- 
Rostoker (KKR) type method developed in Ref. HI for 
the scalar Casimir effect. Because at short distance 
RSPM's provide the dominant contribution to the Cas- 
mir interaction 3 ^' 38 , it would be then necessary to care- 
fully take into account the contributions of all the peri- 
odic orbits associated with the SP. In this context, the 
description that we gave in Sec. IV could be helpful. 
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